Measurement of Lagrangian velocity in fully developed turbulence 
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We have developed a new experimental technique to measure the Lagrangian velocity of tracer 
particles in a turbulent flow, based on ultrasonic Doppler tracking. This method yields a direct 
access to the velocity of a single particle at a turbulent Reynolds number Rx = 740. Its dynamics 
is analyzed with two decades of time resolution, below the Lagrangian correlation time. We observe 



that the Lagrangian velocity spectrum has a Lorentzian form E {u>) 



,Tl/{1 + {TlujY), in 



agreement with a Kolmogorov-like scaling in the inertial range. The probability density function 
(PDF) of the velocity time increments displays a change of shape from quasi- Gaussian a integral 
time scale to stretched exponential tails at the smallest time increments. This intermittency, when 
measured from relative scaling exponents of structure functions, is more pronounced than in the 
Eulerian framework. 



PACS numbers: 47.27.Gs, 43.58.+Z, 02.50.Fz 



Lagrangian characteristics of fluid motion are of fun- 
damental importance in the understanding of transport 
and mixing. It is a natural approach for reacting flows or 
poUutant contamination problems to analyze the motion 
of individual fluid particles j|] . Another characteristic of 
mixing flows is their high degree of turbulence. For prac- 
tical reasons, most of the experimental work concerning 
high Reynolds number flows has been obtained in the 
Eulerian framework. Lagrangian measurements are chal- 
lenging because they involve the tracking of particle tra- 
jectories: enough time resolution, both at small and large 
scales, is required to describe the turbulent fluctuations. 

Early Lagrangian information have been extracted 
from the dispersion of particles, following Taylor's ap- 
proach. Recently numerical and experimental studies 
have focused on resolving the motion of individual fluid 
or tracer particles. The emerging picture is as follows. 
The one-component velocity auto-correlation function is 
quasi-exponential with a characteristic time of the order 
of the energy injection scale ||, Q . The velocity power 
spectrum is expected to have a scaling E^{lo) oc as 
recently reported |j and expected from a Kolmogorov 
similarity arguments. In the same spirit, the second or- 
der structure function should scale as D^^ir) — Coer, 
where e is the the power dissipation. Measurements of 
atmospheric balloons have given Cq = 4 ± 2 and a 
limit Co ^ 7 has been suggested in stochastic models ^ . 
Recent experiments |o| using high speed optical tech- 
niques have shown that the statistics of the Lagrangian 
acceleration are strongly non-Gaussian. 

We have developed a new experimental method, based 
on sonar techniques pd| | , in order to study in a laboratory 
experiment the Lagrangian velocity across the inertial 
range of time scales. We obtain the first measurement of 
single particle velocity for times up to the flow large scale 
turnover time, at high Reynolds number. In this Letter, 



we report the results of this measurements and compare 
with previous observations and numerical predictions. 

Our technique is based on the principle of a continu- 
ous Doppler sonar. A small (2mmx2mnr) emitter con- 
tinuously insonifies the flow with a pure sine wave, at 
frequency /g = 2.5 MHz (in water). The moving particle 
backscatters the ultrasound towards an array of receiv- 
ing transducers, with a Doppler frequency shift related 
to the velocity of the particle: 27rA/ = q.v. The scatter- 
ing wavevector q is equal to the difference between the 
incident and scattered directions. A numerical demod- 
ulation of the time evolution of the Doppler shift gives 
the component of the particle velocity along the scat- 
tering wavevector q. It is performed using a high resolu- 
tion parametric method which relies on an Approximated 
Maximum Likelihood scheme coupled with a generalized 
Kalman filter . The study reported here is made with 
a single array of transducers so that only one Lagrangian 
velocity component is measured. 

The turbulent flow is produced in the gap between two 
counter-rotating discs This setup has the advan- 

tage to generate a strong turbulence in a compact region 
of space, with no mean advection. In this way, parti- 
cles can be tracked during times comparable to the large 
eddy turnover time. Discs of radius R = 9.5 cm are 
used to set water into motion inside a cylindrical vessel 
of height H = 18 cm. To ensure inertial entrainment, 
the discs are fitted with 8 blades with height ht — 5 mm. 
In the measurement reported here, the power input is 
e = 25 W/kg. It is measured on the experiment cool- 
ing system, from the injection-dissipation balance. The 
integral Reynolds number is Re — R^Q,/v — 6.5 10^, 
where VI is the rotation frequency of the discs (7.2 Hz), 
and V = 10~^ m^/s is the kinematic viscosity of wa- 
ter. A conventional turbulent Reynolds number can be 
computed from the measured rms amplitude of veloc- 



2 



ity fluctuations (u^ 



0.98 m/s) and an estimate of 



the Taylor microscale (A = i/lSi^u^^^/e = 0.88 mm); 
we obtain R\ = 740. This value is consistent with earlier 
studies in the same geometry; it corresponds to the range 
of turbulent Reynolds numbers where measurements of 
particle acceleration have been reported . 

The flow is seeded with a small number of neutrally 
buoyant (density 1.06) polystyrene spheres with diame- 
ter d = 250 /im . It is expected that the particles follow 
the fluid motion up to characteristic times of the order 
of the turbulence eddy turnover time, at a scale corre- 
sponding to their diameter, i.e. Tmin d/ud ^ e~^/^d^/'^, 
using standard Kolmogorov phenomenology. For beads 
of diameter 250 /xm, one estimates Tmin ~ 1-3 ms. This 
value is within the resolution of the demodulation algo- 
rithm whose cut-off frequency is at 3 kHz. Note that 
the Kolmogorov dissipative time (t^ — i/i^/e = 0.2 ms) 
is smaller, so that we do not expect to resolve the dis- 
sipative region. The statistical quantities are calculated 
from 3 X 10^ velocity data points, taken at a sampling 
frequency equal to 6500 Hz. The acoustic measurement 
zone is in central region of the flow, 10 cm thick in the 
axial direction and almost spanning the cylinder cross- 
section. In this region the flow is a good approximation 
to isotropic and homogeneous conditions: at all points, 
the mean velocity is non zero, but equal to about one 
tenth of its rms value. 

We first consider the Lagrangian velocity auto- 
correlation function: 



{V{t)v{t + T))t 



(1) 



We observe - Fig. la - that it has a slow decrease which 
can be modeled by an exponential function Pvir) ~ 
e"'^/"^^. This expression defines an integral Lagrangian 
time scale = 22 ms. For comparison, the period of 
rotation of the discs is 140 ms and the sweeping period 
of the blades is 17 ms. The measured Lagrangian time 
scale thus appears as a time characteristic of the energy 
injection. The exponential reproduces extremely well the 
variation of the auto-correlation function, from about 5t,, 
at small scales to AT^. These limits coincide with the 
upper and lower resolution of the technique, so that we 
observe an exponential decay over the entire range of our 
measurement. However, as the variance of the accelera- 
tion must be finite there has to be some lower cut-off 
to this behavior, at times of order r^. These observa- 
tions extend and confirm previous numerical and exper- 
imental studies at moderate Reynolds numbers 1^, H • 
Note that the exponential decay of the Lagrangian veloc- 
ity auto-correlation is a key feature of stochastic models 
of dispersion since it appears as a linear drift term in a 
Langevin model of particle dynamics [|l], ^. 

We show in Fig. lb the velocity power spectrum, com- 
puted both from the data and as the Fourier transform 
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FIG. 1: (a) Velocity autocorrelation function. A best expo- 
nential fit is p^{r) = 1.03e~^^'^'^. It is shown, slightly shifted 
for clarity, as the linear curve in the inset, (b) correspond- 
ing power spectrum; the upper curve is the Lorentzian func- 
tion calculated from the exponential fit of the auto-correlation 
function (shifted for clarity). 



of the exponential decay of the auto-correlation function: 



1 + {Tl^Y 



(2) 



We observe a clear range of power law scaling E'^(uj) cx 
u)^'^. This is in agreement with a Kolmogorov K41 pic- 
ture in which the spectral density at a frequency a; is a 
dimensional function of cu and e: E^[ijj) oc ew~^. To our 
knowledge, this is the first time that it is directly ob- 
served at high Reynolds number and in a laboratory ex- 
periment, although it has been reported in oceanic stud- 
ies 1^ and in lower Reynolds number direct numerical 
simulations |^ . Departure from the Kolmogorov behav- 
ior is observed at low frequencies in agreement with the 
exponential decay of the auto-correlation. At high fre- 
quencies, the spectrum deviates from the Lorentzian form 
due to the particle response. Note in Fig. lb that the 
measurement is made over a dynamical range of about 
60 dB. 

We now consider the second order structure function 
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of the velocity increment 
^2(r)-(Mt + r) 



(3) 



We emphasize that these are time increments, and not 
space increments as in the Eulerian studies. The profile 
D2 (t) is shown in the inset of Fig. 2. It is linked to the 
auto-correlation by I?|'(t) = 2u'^^^ (l — R^{t)) : at small 
times one observes the trivial scaling (t) cx; and at 
large times (r) saturates at 2^^^^^^ (as v{t) and v{t + T) 
are uncorrelated). 
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FIG. 



2: Second order structure function. Inset: profile 
D2{r) as a function of time, non-dimensionalized Tl- In 
the main figure the second order structure function is non- 
dimensionalized by the Kolmogorov scaling er. 

In between these two limits, one expects an inertial 
range of scales with a Kolmogorov-like scaling 



C^er 



(4) 



where Cq is a 'universal' constant. Such a behavior is con- 
sistent with dimensional analysis and with an scal- 
ing range in the velocity power spectrum. Fig. 2 shows 
1^2 (r)/er; a plateau with a constant Co is not observed. 
Note that this also the case in Eulerian measurements 
when the third order structure function is represented in 
linear coordinates . The function reaches a maximum 
at 20t^, for which Co ~ 2.9. This value is in agreement 
with the estimation Cq = 4 ± 2 in and in the range 
of values (between 3 and 7) used in stochastic models for 
particle dispersion [Q. In our case there may also be 
a bias at small times due to particle effects. However if 
we assume the exponential fit for the velocity autocorre- 
lation function to be valid down to the smallest scales, 
we obtain a value Co = 3.5 as an upper bound for the 
maximum of D2 (r)/eT. In our set of measurements be- 
tween R\ = 100 and R\ = 1100, we have observed an 
increase of Co (defined in the same way) from 0.5 to 
4. We point out that in the absence of an equivalent 
of the Karman-Howarth relationship for the Lagrangian 
time increments, a limit value of Co is not a priori fixed. 
Dimensional analysis yields £'2 (''') ^ Co{Re)eT and simi- 
larity arguments give Co{Re) const, or Co(i?e) — * i?e" 
in the limit of infinite Reynolds numbers. 



To further describe the statistics of the Lagrangian ve- 
locity fluctuations, we have analyzed the statistics of the 
velocity increments A^-v. Their PDF 11^ for r covering 
the accessible range of time scales is shown in Fig. 3. 




FIG. 3: PDF (TtIIt of the normalized increment Avt/o-t- 
The curves are shifted for clarity. From top to bottom: 
T = [0.15, 0.3, 0.6, 1.2, 2.5, 5, 10, 20, 40] ms. 

To emphasize the functional form, the velocity incre- 
ments have been normalized by their standard deviation 
so that all PDFs have unit variance. A first observation 
is that the PDFs are symmetric, in agreement with the 
local symmetries this flow. Another is that the PDFs al- 
most Gaussian at integral time scales and progressively 
develop stretched exponential tails for small time incre- 
ments. At the smallest increment, the stretched expo- 
nential shape is in agreement with measurements of the 
PDF of Lagrangian acceleration at identical Reynolds 
numbers ||l^. In our case, the limit form of the velocity 
increments PDF is not as wide as that of the acceleration 
because the Kolmogorov scale is not resolved. Note that 
in regards of the evolution of the PDF, the intermittency 
is at least as developed in the Lagrangian frame as it is 
in the Eulerian one [IT5I . 
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FIG. 4: Evolution of the excess kurtosis factor K{t) = 
((Ar«)'')/((Ar«)^)^ - 3 for the PDFs of the time velocity 
increments. 

The continuous evolution with scale can be quantified 
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using the flatness factor. We show in Fig. 4 the vari- 
ation excess kurtosis K{t) = {{Arv)^) / {(Arv)'^)'^ - 3. 
It is null at integral scale as expected from the Gaus- 
sian shape of the PDF and increases steeply at small 
scales. Below about 5t^, the increase is limited by the 
cut-off of the particle; an extrapolation of the trend to r,, 
yields K{Tjf) ~ 40 in agreement with acceleration mea- 
surements in llOll. 
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FIG. 5: ESS plots of the structure function variation (in dou- 
ble log coordinates). The solid curves are best linear fits with 
slopes equal to = 0.56±0.01, 1.34±0.02, 1.56±0.06, 1.8±0.2 
for p = 1, 3, 4, 5 from top to bottom. Coordinates in arbitrary 
units. 

More generally, one can choose to describe the evolu- 
tion of the PDFs by the behavior of their moments (or 
'structure functions') D^{t) = {\drv\''). Indeed, a con- 
sequence of the change of shape of the PDFs with scale 
is that their moments, as the flatness factor above, vary 
with scale. Classically in the Eulerian picture, one ex- 
pects scaling in the inertial range, (r) oc r^' , at least 
in the limit of very large Reynolds numbers. At the fi- 
nite Reynolds number where most experiments are made, 
the lack of a true inertial range is usually compensated 
by studying the relative scaling of the structure func- 
tions - the ESS ansatz [0. We use the second order 
structure function as a reference. Indeed the dimensional 
estimation of D2 (as that of D^) depends linearly on 
the increment and on the dissipation. Fig. 5 shows that, 
as in the Eulerian frame, a relative scaling is observed 
for the Lagrangian structure functions of orders 1 to 5, 
D^{t) oc I?2 (t)^«. We observe that the relative expo- 
nents follow a sequence close to, but more intermittent 
than the corresponding Eulerian quantity. Indeed, we 
obtain: Cf/Cf = 0.42, ^^/^^ = 0.75, ff/ff = 1.17, 
^5/^3 — 1-28 to be compared to the commonly ac- 



cepted Eulerian values /Cf = 0.36, /Cf = 0.70, 
^i/^i = 1.28, Ci/^i = 1.53. 

In conclusion, using a new experimental technique, 
we have obtained a Lagrangian velocity measurement 
that covers the inertial range of scales. Our results are 
consistent with Kolmogorov-like dimensional predictions 
for second order statistical quantities. At higher orders, 
the observed intermittency is very strong. How the 
Lagrangian intermittency is related to the statistical 
properties of the energy transfers is an open question. 
From a dynamical point of view, the Navicr-Stokes 
equation in Lagrangian coordinates could be modeled 
using stochastic equations. Work is currently underway 
to compare the dynamics of the Lagrangian velocity to 
predictions of Langevin-like models. 
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